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Abstract 

We study the dynamics of ultracold atoms trapped in optical double-well poten- 
tials in presence of noise generated by an external environment. When prepared in a 
Fock number state, the system shows phase coherence in the averaged density profile 
obtained using standard absorption image techniques, that disappears in absence of 
noise. This effect indicates that also in many-body systems an external environment 
may enhance quantum coherence, instead of destroying it. 

1 Introduction 

Cold atoms trapped in a periodic optical potential have become a preferred test ground 
for the study of quantum effects in many-body physics, e.g. in quantum phase transi- 
tion, and matter interference phenomena, and also for applications in quantum infor- 
mation (for instance, see [ll[2] and references therein). 

In particular, the study of phase coherence between fractions of a same condensate 
confined in a double well potential has been the focus of much theoretical and exper- 
imental analysis as a powerful technique to obtain information about the state of a 
trapped condensate. In general, phase coherence cannot be directly observed within 
the trap; nevertheless, it can be put into evidence by releasing the confining potential 
and letting the condensate freely expand, so that the atoms from the two wells may 
interfere with each other: as a consequence, fringe patterns are formed in the spatial 
density. The latter can be obtained by illuminating the expanding clouds with a probe 
light and then by collecting the corresponding absorption image. In particular, infor- 
mation about the initial relative phase can be extracted from the interference pattern 

One expects to observe interference fringes in the density profile only if there is 
phase coherence between the wells. Instead, they appear even starting from a state 
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with random relative phase, hke a Fock number stated In this case however, the in- 
terference patterns have different offset positions in successive shots; as a consequence, 
by averaging over multiple realizations of the experiment the interference pattern dis- 
appears. 

In the following, we shall study how this picture changes in the presence of an ex- 
ternal environment, weakly affecting the dynamics of the condensate in the double- well 
trap. Usually, the environment is a source of noise and dissipation that are expected 
to lead to decoherence. Instead, we shall show that, the environment may create phase 
coherence between two initially uncorrelated fractions of the condensate, an effect that 
may be traced to the presence of an environment induced current between the two 
wells [H]. 

We will consider a system consisting of atoms confined in a double- well trap with 
a very high barrier; in absence of noise, the Fock number states with k atoms in one 
well and N — k in the other are equilibrium states. The presence of an environment 
changes the situation and makes the system an open quantum system; its dynamics is 
no longer unitary and is characterized by the presence of noise and dissipation, so that 
the Fock number states are no longer stationary. We will assume that the interaction 
between the atoms and the environment be weak, with very fast decaying environment 
correlations. In this case, the reduced dynamics of the trapped gas admits a Markovian 
{i.e. memory- less) time-evolution described by a family of maps forming a quantum 
dynamical semigroup \15 \ \16 \ \T7\. An explicit example of such an environment has been 
proposed in [18] , where cold bosons in an optical lattice are coupled to the Bogoliubov 
excitations of another condensate, which mediates transitions between the lowest and 
the first excited lattice Bloch bands. Alternatively, one may consider the coupling of 
the atoms to an external classical stochastic field, with white noise correlations. 

Using open quantum system techniques, we shall analyze below the averaged density 
resulting from the free expansion of the gas after trap release. We shall see that, even 
starting from a Fock number state, with no definite relative phase between the states 
in the two wells, interference fringes may appear in the averaged density profile; they 
disappear in absence of the environment, thus showing the purely noisy origin of the 
effect. We shall also see that this phenomenon may be related to the presence of a 
dissipative current between the two wells [14j . 



2 Cold atoms in a double-well potential 



The usual description of cold atoms trapped in a double well potential uses a two-mode 
Bose-Hubbard model [20\ [2T] with Hamiltonian 

Hbh = £i a{ai + 62 a\a2 + U {a\aif + (0^02)^ - T {a\a2 + aia\) , (1) 



^This phenomenon has been the focus of various investigations, aiming to explain single shot interference 
formation via the detection process (see [6]-[13] and references therein). An alternate explanation is however 
possible: as further discussed in Section 3, it makes use of a projection onto fixed-phase, coherent-like states. 
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where ai 2, a| 2 annihilate and create atom states in the first, second well, respectively, 

and satisfy the Bose commutation relations [oj , aj] = 5ij. Of the contributions to Hbh, 
the last one corresponds to a hopping term depending on the tunneling amplitude T, 
the first two are due to the trapping potential and are proportional to the depth of 
the wells: in the following we will consider the case of a symmetric trap, i.e. £1 = £2- 
Finally, the third term, quadratic in the number operators a|aj, takes into account 
repulsive Coulomb interactions inside each well. 

Let {|<^i)}^i be a complete set of orthonormal single-particle atom states; the 
creation operator ip^{x) of an atom at position x can then be decomposed as 

00 

^t(x) = ^<^:(x)al , (2) 

i=l 



where a| creates an atom in the state |</)i) = a||uac) with wavefunction (j)i{x) 
The Bose-Hubbard Hamiltonian ([1]) results from a tight binding approximation, where 
only the first two of the basis vector are relevant; in this case (pi^2ix) are orthogonal 
functions, localized within the first well, (/>2 within the second one. 

The total number of atoms is conserved by ([T]). Therefore, the Hilbert space of 
the system is + 1-dimensional and can be spanned by Fock number states of the 
form 

\k,N-k) := K)' ("2r~' l^ac) , (3) 



/fc!^(Af-A;)! 

describing the situation in which the first well is filled with k atoms, while the other 
one contains N — k particles. They are obtained by the action of the creation operators 
on the vacuum state. These states turn out to be eigenstates of the Bose-Hubbard 
Hamiltonian when the tunneling term can be neglected. 
Alternatively, one can consider the so-called phase-states 



\V,C;N) = [y^e'^al + yr^e-^^aj)"^ \vac) , (4) 



depending on two real parameters, an amplitude ^ G [0, 1] and a phase if G [0, 27r]. 
These states describe a physical situation where all A^ atoms are in a same superposition 
of a two-mode single particle state with relative phase ip ; from a many-body point 
of view, such states correspond to a macroscopic wavefunction obeying the Gross- 
Pitaevski equation [3]. Moreover, the parameter ^ is related to the mean value of the 
relative occupation number, (99, ^; A^| (ajai — 02a2)|V;C;^) = -^(2^ — !)• 
These states form an overcomplete set 

/\r _|_ 1 fi /■27r 
1^ = ^- de d^\^,i;N){ip,i-N\. (5) 

Although not orthogonal, for sufficiently large A^ (a very common situation in actual 
experiments), the phase states become nearly so, 

{if, C; N\^',e;N)^^ - e) - . (6) 
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One can then decompose any N atom state in terms of these states; in particular, for 
Fock number states, one can write 

/V _|_ 1 /■! /"St 

\k,N-k) = ^— d^J^ d^{^,^;N\k,N-k)\^,C;N) , (7) 
where the overlapping functions are explicitly given by 

^; N\k,N-k)= (1 - e-^e^-f ) . (8) 

Further, the action of ip{x) on the phase states can be easily derived using 

3 Density profile after free expansion 

As explained in the Introduction, our aim is to study the spatial density profile obtained 
by releasing the trap and by subsequently letting the cold atoms freely expand for a 
time T. If r is large enough, the two fractions of cold atoms interfere with each other 
thus providing a spatial fringe-pattern that one can probe by light-absorption. 

Suppose the trap is released at r = when the cold atoms are in the phase state 
In order to describe the free expansion, disregarding atom-atom interactions, we 
make use of the following model: we let the single-particle states evolve under the 
action of an unitary operator Ur', notice that, by the unitarity of Ur, the basis states 
|(/>j) evolve in another set of basis states, 

\<Pi), = Ur\(|>^) := al(T)\vac) , aj(r) = Ur^Ul . (10) 

Thanks to the absence of interactions during the expansion, the evolution of the 
phase states follows immediately from the evolution of the single particle states, 

|<^,C;iV>^ = (x/^e'^al(T) + yr^e-^^4(r))'^ \vac) , (11) 

while both the relations ([5]) and 1^ still hold for any r [22] . 

The annihilation operator (which does not evolve in time) could be again decom- 
posed over the new basis 

oo 

1p{x) ='^(pi{x,T)ai{T) , (pi{x,T) = {x\(j)i)r , (12) 

i=0 

and therefore the action on the freely evolved phase state reads 

i^ix) I ^, e; N)^ = Vn{^<P, (x, t) e'^ + y^T^ r) e"^^) \(l),i;N - 1>^ . (13) 

If the two fractions of the condensate expand independently, we can explicitly 
evaluate the relative phase gained during the ballistic expansion. Writing (pi{x,T) = 
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r)|e*^*(^''^), the phase 9i{x, r) can be obtained, for large r, in terms of the classical 

behavior of the velocity v = — V^j = — [2j- We shall focus upon the direction x, hence 

m T 

upon the phase modulation 

lm(x-xi,2)^ 

6'i,2(a;,r)- 



2 hT 

where Xi is the position of the i-th well; it follows that the overall relative phase after 
a time of flight r can be well approximated by 

Ttl u 

0,{x,r)-e2{x,T) = —x, (14) 

hT 

where the trap is considered to be located along the x axis, centered at x = with the 
two wells at position xi = —d/2, X2 = d/2, respectively. 

Then, starting from a phase state, the one-particle spatial density explicitly reads 

= r{v,i;N\'4j\x)ij{x)\ip,i;N)^ 

= iv(^el0i(x,r)|2 + (l-0|</'2(x,r)|2 + 

/ md 



+ 2Ve(l-e)|</)i(:E,T)| |02(x,t)| cos ^—x + v9j j . (15) 

The last term in this equation exhibits the presence of fringes that are equally separated 
hr 

by a distance (. oc — - with an offset position fixed by the initial phase ip. 
md 

The situation is completely different if one starts instead from a Fock number 
state; in this case the mean value of the density operator after the free expansion can 
be calculated by decomposing the state using the analog of ([7]) at time r, 

\k,N-k)^ = ^— dU d^{^,^;N\k,N -k)\^,C;N)^ . (16) 
^"^ Jo Jo 

Experimentally, the light-absorption image relative to a Fock number state shows 
the same fringe pattern (jl5p of the superfluid state, whereas one would expect no 
interference fringes at all. However, if one repeats the experiment with the same initial 
conditions, the fringe patterns are obtained with randomly distributed offset positions. 
Therefore, by averaging over them, the fringes cancel out and the averaged density 
profile does not exhibit any interference pattern [3]-|13j. 

In order to interpret this effect, notice that the absorption image amounts to a 
measurement of the initial relative phase ip, for it gives a one-shot density profile of the 
form ()15|) . Furthermore, starting from a phase-state ^; A^), every single repetition 
of the experiment will always give the same interference pattern n^^g(x,r) [11]. On 
the contrary, when the initial state at r = is a Fock number stare \k, N — k), each 



^Note that tlie scalar products {ip, ^; N\k, N — k) in pB]) do not depend on time. 
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single shot will correspond to a random selection of one of the phase-states (jlip in the 
expansion (|16p . with probability densities that are independent of the relative phase 
and proportional to 

i((^,e;iv|A;,iv-fc)|2= (^^^e'(i-e)^-' • (17) 

In other words, every single-shot measurement of the density profile in the state 
|/c, — A;)^ resulting from the Fock number state \k,N—k^, randomly projects out the 
density profile of one of the phase-states jv?, ^; -^)^, with weight (fT7|) . As a consequence, 
the averaged density profile reads 

N + 1 f^'^ 
nk{x,T) := — / dU dif\{^,^;N\k,N -k)\'^n^^i:{x,T) 

Jo Jo 

^ \k + l)\Mx,r)\^ + {N -k + l)\M^,T)\^] , (18) 



iV + 2 



and shows no interference pattern (see |22| for further details). 

The density profile at time r after trap release depends on the state of the A^-atom 
system at time r = 0. Since we are interested in studying the effects of a dissipative 
open dynamics during the interval of time t spent by the system within the confining 
double-well potential, we will describe the system state at time t by a generic density 
matrix pj: 

N 

Pt=Yl RUk,N-k){q,N-q\ . (19) 

k,q=0 

Notice that the information about the dynamics during the time t spent by the atoms 
within the double- well trap is contained in the coefficients i?^^. 

In order to construct the averaged density profile in such a case, one considers 
the ensemble of profiles n^^^{x,T) = T-(^(p, S,; N\'iIj^ {x)'i(j{x)\(p, S^; N'j that are obtained 
in many repetitions of the experiment, with the same initial conditions. Each image 
shows interference fringes of the form (jlSp and appears in the statistical ensemble of 
all collected images with weights proportional to 



{^,^;N\p,\^,^;N) = f; i?^,./(f)ne^(l-e)^-^e'-(^-) . (20) 
k,q=o V ^ / V y / 

The resulting one-particle spatial density profile is then given by the corresponding 
statistical average 



N + 1 







2tt 



N 



Y,Rik {k + i)\Mx,T)\^ + {N-k + i)\Mx,T)\' 



k=0 



N + 2 

N-l ^ 

+ ^Wi + l)(iV - A;) |0i(x,r)| \Mx,t)\ e"^^^ + c.c. . (21) 



k=0 
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The presence of visible interference fringes is strictly related to the next-to-diagonal 
entries of the density matrix Rl.f^_^_i (and their hermitian conjugates) H 

Notice that the previous average can be computed as the mean value of the density 
operator rl^^{x)ip{x) with respect to a density matrix pt which results from the freely 
evolved state Ur pt ut by the action of a generalized measurement, a so-called Posi- 
tive Operator Valued Measure (POVM) [23\ I24j. consisting of the family of projectors 
P<fi,C,N ■= \^,(,] ]^){^, S,; N\. Explicitly, one has 



N + 1 
2tt 



rl 1-211 

d(f P^,^;N Ur Pt Ul P^,^-N ■ 

Jo Jo 



(22) 



The map from pt to pt corresponds to the canonical description of the quantum mea- 
surement process associated with the POVM {P^^^.jv}- One notices that it does not 
correspond to a measurement of the number operator in the state Ur pt Ul which gives 
an average number of the form 



TV 



)-0(j;) UrPtUfj = N dx2- ■■xn{x, X2, . . . , XN\UrPtUl\x, X2, . . . , Xn) 



N 



k=0 



N-l 



+ ^{k + 1){N- k) \<i,i{x, t)\ \ct>2{^, t)\ e- 



+ c.c. \ , (23) 



A:=0 

which is often used to fit experimental data; this average corresponds to a POVM 
consisting of projections onto eigenstates of the operator 'ip'^ {x)iIj{x), not onto phase- 
states. 

For large, the difference between (j23p and (|2ip results very small for all practical 
purposes. Furthermore, the difference does not show up in the oscillating terms. In 
the following we will make use of the the description of the absorption images in terms 
of a POVM consisting of the projections onto the phase states, since this appears to 
describe more adequately the absorption image formation 



4 Noisy dynamics within the trap 

In this Section we study the dynamics of the atoms within a trap immersed in an 
external environment; in particular, we shall focus upon its effects on the ground state 
|A^/2, A^/2) of ([1]) with a very high potential barrier. 

The density matrix of the total system, atoms plus environment, psE evolves uni- 
tarily under the action of a total hamiltonian which can be written in the form 

HsE = Hbh ^1e + I^He + Hi , (24) 

^ The evolution of the atoms within the trap affects only the coefficient Rj,^ and depends on the time t 
spent inside the trap; the time of flight r refers instead to the time interval between the releasing of the trap 
and the snapshot of the density and it is a parameter of the experimental setup. 
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where H^h is as in ([T]) and describes the motion of the system in absence of the envi- 
ronment, He is the hamiltonian of the environment, while Hj describes the interaction 
between them and can be taken to have the general form 

Hi = Y,yi®Bi (25) 

i 

where Vi and Bi are suitable system and environment hermitian operators, respectively. 

The evolution of the system density matrix p = Ti e{pse) can be obtained by 
tracing out the environment degrees of freedom from the standard unitary evolution of 
the total system initial state psE- In general this operation leads to a master equation 
affected by memory terms and nonlinearities. The latter can be avoided by choosing 
an initially uncorrelated state of the system and the bath, psE = P ® Pe] this is a 
common situation in experimental contexts when the environment is supposed to be 
in a reference equilibrium state pE ■ 

On the other hand, memory effects can be considered a short time phenomenon; 
if we are not interested in this transient regime, they can be neglected by studying 
the effective evolution of the system on a slower timescale. This amounts to a Markov 
approximation |161ll5i[T7] which is only plausible when a neat separation exists between 
system and environment timescales. Necessary conditions for this approximation are 
for instance satisfied within the so called singular coupling limit |25[ [26l [27] , where the 
characteristic decay times of the environment correlations, 

Gi,{t) = {Bi{t)Bj) = Tr [Bi{t) B, pe] , Bi{t) = e''"^Bi e"'*^^ , (26) 

are very small when compared to the typical time scale of the system. This condition 
is satisfied by a thermal bath in equilibrium at very high temperature or by an external 
stochastic white-noise classical field; in this latter case, the matrix [Gij{t)] turn out 
to be real symmetric. In [18], a suitable bath fulfilling the above conditions has been 
proposed, where the system is coupled with the Bogoliubov excitations of another 
condensate, which drive a transition of the trapped atoms to the first excited lattice 
Bloch band, with a characteristic transition rate k0 

In the singular-coupling limit, the dissipative time-evolution of the open system 
density matrix is generated by a master equation of the form 

Ml = L[p(t)] ^ -z[Hbh + H^'\ Pit)] + B[pit)] , (27) 

where the noise contributes with an effective correction H^"^^ to the free hamiltonian 
and with a linear operator B[-] which cannot be recast in hamiltonian form. Their 



is interesting to notice that another possible mechanism inducing dissipation is provided by inelastic 
collisions among atoms on the same lattice site, that can be tuned on by a nonvanishing imaginay part of the 
interaction parameter U [19]. These effects would nevertheless lead to atom loss, i.e. to non conservation of 
TV, requiring an extension of the formalism to be properly treated. 
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explicit expressions take the standard Kossakowski-Lindblad form [151 HSl IZI : 



ij=l 
ij=l 



(28) 



(29) 



The operators Vi are the operators involved in the interaction hamiltonian (j25p . while 
the entries of the two hermitian matrices [cij] and [sij] embody all the information 
about the effects induced by the environment and are obtained from the correlation 
functions of the environment: 

1 







dt Gij{t) 



27r 



uj=0 



dio ■ 



UJ 



(30) 
(31) 



where J- denotes Fourier transforming and V the principal value. 

In particular Cjj is known as Kossakowski matrix and turns out to be positive 
definite, a key requirement to ensure the positivity of the eigenvalues of the evolved 
density matrix pt at all times and therefore the consistency of the dynamics. For 
instance, in the case of the environment proposed in [18], the matrix [cjj] takes the 
explicit expressiorll 

-2k 
2k 




[Cij] 



( 2k 
-2k 


-in 










IK \ 

— IK 



k/2/ 



(32) 



The operator L[-] in ()27p generates a family of maps 7^ = e , which forms a 
quantum dynamical semigroup \15 \ I16[ [T7| . fulfilling the forward in time composition 
law jt ° 7s = 7t+s) t,s > 0, and describing the irreversible character of the dissipative 
evolution. 

Quite in general, the presence of an environment will affect both the tunneling 
amplitude and the minima of the confining potential; the environment operators Bi 
will then be coupled to the following bilinear atom operators appearing in the system 
hamiltonian ([T]) a\ai, 0302, a\a2 and CgOi, or more precisely to the four hermitian 
operators 

Vi = , a\a2 , {a\a2 + a\ai) , i {a\a2 — a\ai)^ . (33) 

As a final remark, notice that the action of the environment ceases as a consequence 
of the trap release. 



^The dissipative effects induced by inelastic atom collision discussed in [19] can similarly be described in 
terms of a master equation of the form (|27p : the magnitude of the elements of the corresponding Kossakowski 
matrix is determined by Tm{U). 
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5 Effects of noise in the density profile 



As mentioned before, we shall consider an initial atom state of the form 

po = \N/2,N/2){N/2,N/2\. (34) 

In absence of the environment and for a high potential barrier, this state is the ground 
state of ([I]); as a consequence, the components R\q in (fT9l) are time independent and 
Rkq = if A: and q differ from N/2. Therefore in both (|2ip and (|23p no interference 
fringes appear. 

In order to study the effects induced by the presence of noise, it is sufficient to 
consider the small time dynamics resulting from (j27p . whose solution up to first order 
in time then reads 



Pt = Po 



+ it [Hbh + H^^\ Po] +mpo] + 0{t^) . (35) 



Since [Hbh, Po] = 0, only the dissipative contributions in ([35]) need to be studied. 
The explicit expressions of ©[/jq] and i[H^'^\pQ] are collected in the Appendix. Using 
them, to ffi'st order in t, one finally gets 

npAx,T) = y (|0i(x,r)|2 + |02(x,T)|2) + 

+ t — <^4 9m(c34)(|02(x,r)|2-|0i(x,r)|2) + 



4 

+ \Mx,t)\\Mx,t)\ 



^2 3f?e(s24 — sii) + 9m(ci4 — 024)^ cos 



md 



+ (2 3f?e(s23 - S13) + 9m(ci3 - 023)^ sin | - — x 



/md 



(36) 



The last two lines of the last equation contain oscillating terms which give rise to 

md ' 



interference fringes spaced by a distance £ oc notice that these fringes disappear 



in absence of noise {i.e. when Cij = = Sij). While the interference pattern is the 
same as in ()15p and is essentially due to the ballistic expansion, the amplitude of these 
oscillations is proportional to the entries of the matrices [cij] and [sjj], which depend 
on the strength of the environment correlation functions. 

Since for large enough r the atomic clouds coming from the two wells overlap, i.e. 
\(Pi{x,t)\ ~ \4>2{x,t)\, the magnitude of the dissipative terms goes as tN\cij\ with 
respect to the standard, noise independent contribution. Despite the large number 
of atoms, because the system environment coupling is weak and thus the constants Cij 
and Sij are small [25] [26l [27] , the first order expansion in (|35p is meaningful even for 
sufficiently long system environment interaction times t. 

For instance, if the Kossakowski matrix [cij] is as in (|32p . the magnitude of its 
entries is proportional to the transition rate k; for typical values of this parameter 
and of N, t can be large enough to allow a direct experimental observation of the 
interference fringes. Finally, notice that if instead the system is coupled to a stochastic 
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classical field, the matrix [cj-,] is real symmetric while all the entries of [sij] vanishes [25]; 
therefore, the oscillating terms in (|36p disappear; this means that classical correlations 
do not lead to any environment induced interference phenomena. 



6 Relations with dissipative current 



To better understand the physical origin of the environment induced interference pat- 
tern, let us consider the mean value of the following current operator 



J = i{a[a2 — a\ai) , 



(37) 



which describes the motion of the barycenter of the atom gas in the double well. Its 
mean value is clearly zero in the initial state po, while, in the evolved state up to first 
order in t, it becomes 



(J) 



pt 



■N 
iVl- + l 



^(c23 - cis) + 2 3f?e(si3 - S23) 



(38) 



By comparing the last expression with ()36p one notices that, if the current operator has 
nonzero mean value, then interference fringes can be observed; however the structure of 
the interference pattern is in general richer, as also the tunneling operator a\a2 + a^ai 
contributes to it. Indeed, one finds 



{a\a2 + alai] 



Pt 



■ N 



(ci4 - C24) + 2 3f?e(s24 - S14) 



(39) 



This is precisely the situation that occurs for the Kossakowski matrix (132^ : it leads to 
a vanishing mean value of the current operator, while instead the mean value of the 
tunneling operator is proportional to k. 



7 Outlook 

We have seen that, even in the ground state pq of the hamiltonian ([T]), the presence of a 
noisy environment may give rise to interference effects in the averaged density profile. 
The periodic structure of these fringes comes from the ballistic evolution of the gas 
and does not depend on the noise; their oscillations amplitude are, however, directly 
proportional to the characteristic parameters of the noise and vanish in absence of it. 

This phenomenon can be related to a dissipative current between the wells, so that 
the study of the interference fringes in the averaged density profile gives information 
about this current. 

The system we analyzed can in principle be engineered in the laboratory, for in- 
stance using a second condensed gas as an environment as proposed in [18j : the result- 
ing interferences fringes can be detected by subtracting the density pattern in absence 
of noise. This effect can also be useful to further understand the role of noise in double 
well potentials and in particular for the possibility of driving the system out of the 
ground state po without lowering the barrier. 
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Similar conclusions can be easily extended to the case of an optical lattice; in 
this situation the capability of suitably engineered environment to create correlations 
between different sites of the lattice could be practically very useful in order to create 
specific classes of correlated states. 



Appendix 

By setting |A;) = \k,N — k) for sake of simplicity, the explicit expressions of ©[po] and 
i[H^'^\po] read 

= TIT + ^J (C33 + C44 + ^(C34-Cyj + + 1| 

/ , * AiN .,N . , , |iV, ,iV, 

+ (^C33 + C44 - ^ (C34 - C34) J ly - l)(y - 1| - 2 (C33 + C44) I y/Vy I 



+ < 



TV 1 

~r (^13 — + C23 — i C24) — - {c{-i ~ ^ — C23 + i C24) 



- ^ (ci3 -ICU + C23 - ^ C24 J 



I 2 X 2+^1 



In (N T 



^ (Ci3 + i Ci4 + C23 + i C24) + ^ (Ci3 + i Ci4 - C23 - i C24) 



- (Cl3 + t Ci4 + C23 + ^ C24 J 



I 2 2-^1 



2X2+^1 + 1 2X2- 'I 



+y (y + 1) (c33 - C44 + i (C34 + c*u)) I y + l><y - 1| + /i. c. I 



(40) 
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